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ABSTRACT 

We present an analytical description of the energetics of the population of cosmic accretion shocks, 
for a concordance cosmology (Q m +SI& = 1). We calculate how the shock-processed accretion power 
and mass current are distributed among different shock Mach numbers, and how they evolve with 
cosmic time. We calculate the cumulative energy input of cosmic accretion shocks of any Mach 
number to the intergalactic medium as a function of redshift, and we compare it with the energy 
output of supernova explosions as well as with the energy input required to reionize the universe. 
In addition, we investigate and quantify the effect of environmental factors, such as local clustering 
properties and filament preheating on the statistical properties of these shocks. We find that the 
energy processed by accretion shocks is higher than the supernova energy output for all z < 3 and 
that it becomes more than an order of magnitude higher in the local universe. The energy processed 
by accretion shocks alone becomes comparable to the energy required to reionize the universe by 
z ~ 3.5. Finally, we establish both qualitative and quantitatively that both local clustering as well as 
filament compression and preheating are important factors in determining the statistical properties of 
the cosmic accretion shock population. 

Subject headings: shock waves - large-scale structure of universe - galaxies: clusters: general - inter- 
galactic medium - methods: analytical 



1. INTRODUCTION 

The formation of shocks in the baryonic component 
of matter in the universe is an inevitable and integral 
part of the process of cosmological structure formation. 
The processing of gas through cosmic shocks plays a ma- 
jor role in the distribution of cosmic baryons and their 
energetics and thermodynamic properties. In particu- 
lar, the initial shock heating of primordial gas leaves 
it optically invisible and thus contributes to the inven- 
tory of dark baryons (Hegyi & Olive 1986; Fukugita et 
al. 1998; Spergel et al. 2003; Cyburt et al. 2003). A 
full observational accounting of dark baryons in all of 
their states is now emerging, with the first evidence for 
a large component of diffuse, low-density intergalactic 
gas distributed in the filaments and sheets that com- 
prise the "cosmic web," and detected in the X-ray for- 
est via absorption lines from highly-ionized metals (Hell- 
sten et al. 1998; Perna & Loeb 1998; Fang et al. 2002; 
Nicastro et al. 2002; Finoguenov et al. 2003; Mathur et 
al. 2003; Nicastro et al. 2005). This warm-hot intergalac- 
tic medium is thought to arise in structure formation 
shocks (e.g. Hellsten et al. 1998; Cen & Ostriker 1999; 
Dave et al. 2001; Nath & Silk 2001; Furlanetto & Loeb 
2004; Kang et al. 2005). 

Moreover, structure formation shocks heat the inter- 
galactic medium and are likely to act as acceleration 
sites for nonthermal particles, the "structure formation 
cosmic rays" (e.g. Miniati et al. 2001a, b; Brunetti et 
al. 2001; Fujita & Sarazin 2001; Berrington & Dermer 
2003; Gabici & Blasi 2003a; Brunetti et al. 2004; Kang 
& Jones 2005 and references therein). Recent detec- 
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tions of nonthermal emission from galaxy clusters (e.g. 
Fusco-Femiano et al. 1998; Govoni et al. 2001; Bagchi 
et al. 2002; Feretti et al. 2004 and references therein) is 
consistent with inverse Compton emission by relativistic 
intracluster electrons. Such a cosmic ray population pre- 
sumably also includes hadrons, and so would also have 
distinct 7-ray and light-element signatures, which are 
currently subjects of intense investigation (Loeb & Wax- 
man 2000; Totani & Kitayama 2000; Miniati 2002; Scharf 
& Mukherjee 2002; Suzuki & Inoue 2002; Totani & Inoue 
2002; Berrington & Dermer 2003; Gabici & Blasi 2003b, 
Keshet et al. 2003; Miniati 2003; Reimer et al. 2003; 
Gabici & Blasi 2004; Prodanovic & Fields 2004; Pro- 
danovic & Fields 2005; Suzuki & Inoue 2004; Kuo et 
al. 2005). 

Cosmic shocks occur during different facets of struc- 
ture formation and in a variety of environments, hence 
there are at least two distinct ways to categorize them: 
according to the physical processes driving them, and ac- 
cording to the state of the medium in which they form. 

There are three principal processes associated with cos- 
mic structure formation which result to the formation of 
large-scale shocks. 

1. Accretion of intergalactic gas by a collapsed, viri- 
alized structure. In this case, an accretion shock 
is formed at the interface between virialized and 
diffuse gas (e.g. Bcrtschinger 1985a; Miniati et 
al. 2000) The shock is driven by the gravitational 
attraction exerted on the diffuse gas by the accretor 
and heats the accreted gas to the virial tempera- 
ture. 

2. Merger of two collapsed structures. In this case, 
a merger shock is formed at the interface between 
the gas components of the merging objects (e.g. 
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Miniati et al. 2000; Gabici & Blasi 2003). The 
shock is driven by the mutual gravitational attrac- 
tion between the objects. 

3. Accumulation of intergalactic gas onto large-scale 
cosmic filaments and sheets. Contrary to virial- 
ized structures, large-scale filaments and sheets, al- 
though overdense, are still expanding. At the edges 
between filaments and their surroundings, shocks 
form, driven by the difference in expansion velocity 
between the overdense filament and the underdense 
neighboring voids (e.g. Bertschinger 2005b). More 
complicated shock surfaces inside the filaments are 
additionally produced due to the gravity of the col- 
lapsed objects residing in the filament interior. 

Shocks can also be divided according to the state of 
the gas passing through them, into external and internal 
shocks (e.g. Miniati et al. 2000; Ryu et al. 2003). Ex- 
ternal shocks process pristine material, which has never 
been shocked before by any of the processes described 
above. External shocks are mostly filament shocks, since 
the process of formation of individual virialized struc- 
tures (associated with the other two types of shocks) oc- 
curs principally within filaments, and therefore in most 
cases involves gas which has already been processed at 
least by filament shocks. Internal shocks process gas 
which has already been shock- heated in the past. All 
merger shocks, as well as many accretion shocks, are in- 
ternal shocks. 

Since external shocks process colder material of lower 
sound speed, their Mach numbers are generally higher 
than those of internal shocks. However, because the gas 
passing through internal shocks has already been com- 
pressed, internal shocks process more mass and kinetic 
energy than external shocks. 

In this work, we present an analytical study of the ener- 
getics and cosmic evolution of the population of of cosmic 
accretion shocks, in a concordance f2 m + f2A = 1 universe. 
In addition, we assess the effect of environmental factors 
on the properties of these shocks. The environmental 
factors we consider are: (a) variations in the local mat- 
ter density and temperature of the region in which each 
accretor is embedded imprinted in the primordial density 
field and (b) filament preheating and compression. 

In order to distinguish between the effects of each en- 
vironmental factor, we start from a base model which 
accounts for no environmental effects and regards all ob- 
ject as being embedded in an environment well repre- 
sented by the background universe. In this case, shock 
properties depend solely on the mass of the accretor. To 
this model we then add each of the environmental factors 
mentioned above, each version of our model being refined 
by one effect with respect to the previous version. The 
final (full) model we construct in this way accounts for 
the effects of the mass distribution of accretors as well as 
both environmental factors considered here. 

Although each effect is treated under a number of sim- 
plifying assumptions which necessarily introduce some 
error in our calculations, we make sure to on one hand 
quantify this error whenever possible and on the other 
hand to consistently err on the side of underestimating 
the effect of environment. In this sense, our results are an 
absolute lower limit of the effect of environmental factors 
on the properties of the population of cosmic accretion 



shocks. In addition, the trends and qualitative changes 
introduced to the properties of the accretion shock pop- 
ulation are clearly identified. Any future quantitative 
correction to our model will only enhance such trends. 

The shock properties as a function of object mass are 
derived from the temperature jump across the shock sur- 
face (as is the case in cosmological numerical simula- 
tions), while the underlying mass-function of the accre- 
tors is taken to be the Press-Schechter mass function. 
We calculate statistical quantities characteristic of the 
population of cosmic accretion shocks, such as the dis- 
tribution of number density of objects, processed kinetic 
energy and shocked mass with respect to Mach num- 
ber, as well as their evolution with cosmic time, and we 
use our findings to discuss the relative importance of cos- 
mic accretion shocks in the energetics of the intergalactic 
medium. 

Our paper is organized as follows. In Sj2]we present the 
formalism describing the properties of a single accretion 
shock around a cosmic structure. In fj3]wc combine our 
single-shock model with different cosmological distribu- 
tion functions to derive the statistical properties of the 
population of cosmic accretion shocks. In Sj?]we present 
the results of our model for a concordance f2 m + S1a = I 
universe. In Sj5]we discuss the differences and similarities 
of our approach from other analytic cosmic shock mod- 
els based on the Press-Schechter approach. Finally, we 
summarize and discuss our findings in <JS] 

2. PROPERTIES OF A SINGLE SHOCK 

In this work, the properties of the accretion shock 
around a single virialized host object are derived from the 
temperature jump across the shock, in agreement with 
the method used to derive shock properties in cosmo- 
logical simulations. The pre-shock temperature is taken 
to be an appropriately defined environmental tempera- 
ture, while the post-shock temperature is taken to be the 
virial temperature of the host structure. However, the 
virial temperature is a mean quantity of the host object 
and does not necessarily characterize the temperature 
right behind the shock. To assess the error introduced in 
our calculation of the Mach number due to this assump- 
tion, we compared our temperature-jump technique with 
the Bertschinger (1985) similarity solution which does 
take into account the radial dependences of the post- 
shock gas. To ensure an appropriate comparison, we 
applied our formalism to an Einstein-de Sitter universe 
and in the high-Mach limit, where the Bertschinger solu- 
tion is applicable. We found that the deviation from the 
Bertschinger solution is at the level of 1%. 

Throughout this paper, we assume an adiabatic equa- 
tion of state, and we consider all shocks to be non- 
radiative. We also assume that any individual collapsed 
object as well as its accretion shock are spherically sym- 
metric. We take the accretion shock position around each 
structure to coincide with the virial radius of each struc- 
ture. 

The Mach number of a shock, A4, is defined as the ratio 
of the velocity of the accreted material in the shock frame 
to the adiabatic sound speed of the accreted material. 
The Mach number is related to the temperature jump 
across the shock through 

T 2 _ (5X 2 -l)(A4 2 + 3) 

T x 16M 2 ■ [ ' 
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where T\ and T2 are the pre-shock and post-shock tem- 
peratures correspondingly, and we have assumed a ratio 
of specific heats 7 = 5/3, constant across the shock. In 
the limit Ai ^s> 1 this equation is further simplified, Ai — 
J (I6T2) / (5Ti). The pre-shock temperature can be writ- 
ten in terms of the adiabatic sound speed of the pre-shock 
material c s i as kT\ = /iTOpC^/7, where k is the Boltz- 
mann constant, p, is the mean molecular weight of the ac- 
creted gas, and m p is the proton mass. If we also take T2 
to be the virial temperature of the accretor which has a 
mass m, T vir = ^TO p Gm 2 / 3 (47r/ c p m! o) 1 / 3 (l-|-z)/(5fc3 1 / 3 ), 
then the ratio T2/T1 becomes 



T 2 (m, z) 



2.7x10^ 



fc 



18tt 2 mi 



1/3 

(1+*) 



15 kms 1 



where f c is the compression factor for a virialized ob- 
ject (which may vary with virialization redshift, depend- 
ing on the cosmological model), z is the virialization 
redshift, h is the dimensionless Hubble parameter, and 
mg = 5.96 x lO 14 /i _1 il m M0 is the mass included in a 
sphere of comoving radius r% = 8/i _1 Mpc assuming the 
mean matter density inside the sphere to be equal to the 
cosmic mean. 

The (baryonic) mass current, defined as the rate at 
which mass crosses the surface of a single accretion shock 
around a structure of mass m at an epoch z, is 

Ji = ^ L ^=ATrr 2 v (m)n M (l + z) 3 (l + S s )Mc sl (3) 

where r v is the virial radius of the structure, 

1/3 / , \ -1/3 



fc 



(1 + *) , (4) 



r v = lAh- 1 Mpc — 

p c fi is the critical density at the present cosmic epoch, 
^bA=,o(l + z ) 3 is the cosmic baryon density at the epoch 
of interest, and (1 + J S ) is the density enhancement (with 
respect to the cosmic mean) just outside the shock. 

Finally, the kinetic power crossing the accretion shock 
around a single structure of mass m is P\ = dE/dt = 
0.5(dm/dt)vf = 0.5 Jiw 2 , or 



Pi = 2vrr v (m) 2 O b p c ,o(l + z) 3 (l + S a )M 3 c 3 sl 

3. PROPERTIES OF THE POPULATION OF COSMIC 
ACCRETION SHOCKS 



(5) 



We will examine three distinct models for the cosmic 
accretion shock population. We will start from a "base" 
no-environmental-effects model, in which the properties 
of each accretion shock are simply a function of the accre- 
tor mass. We will then gradually build up more complex 
models by adding one-by-one the different environmental 
factors we wish to examine (local density and tempera- 
ture variations originating in the primordial density field, 
and filament compression and preheating). In this way, 
we will assess the effect of each factor on the properties 
of the cosmic shock population. Our third, most com- 
plex model, will include the effect of both environmental 
factors we are addressing here. 

3.1. Model 1 (base model): no environmental effects 

If we assume that all collapsed objects accrete baryons 
of uniform density and temperature (the density and 



temperature corresponding to the mean, background- 
universe values of diffuse baryons), then at a given red- 
shift, all objects of a certain mass will have accretion 
shocks with identical Mach numbers, and will process 
the same amount of mass and kinetic power per unit 
time. In the relations of £j2l <5 S = and c s i = 15 km/s as 
we will be treating only post-reionization redshifts. The 
mass distribution of collapsed objects will be the one de- 
scribed by the Press-Schechter mass function (Press & 
Schechter 1974; Lacey & Cole 1993) 

1 2 



dn 
dm 



(m, z) — 



2 p. 



m,0 



diner 



dint 



exp 



S c (z) 



2[a(m)Y 



(6) 

(2) where S c (z) is the linearly extrapolated overdensity of an 



object which collapses at redshift z, er(m) = \J S(m) is 
the square root of the variance of the linearly extrapo- 
lated field smoothed at a mass scale m, and p m ,o is the 
cosmic mean matter density at the present time. 

3.2. Model 2: model 1 + effect of primordial density 

field 

In the second variation of our model, we will relax the 
assumption that the temperature of the accreted material 
and the pre-shock density is the same for all structures, 
and we will consider the effect of variations in the local 
environment of accretors of the same mass caused by the 
structure of the primordial density field. In this case, 
accretors of the same mass will reside in a distribution 
of environments, parametrized by the local overdensity 
or undcrdensity (as determined by the evolution of pri- 
mordial density fluctuations), 5$ = (pi OC ai — Pm)/Pm- In 
addition to affecting the pre-shock density enhancement, 
5 S + 1 = Sg + 1, such density variations will also cause 
temperature variations of the pre-shock gas, due to adia- 
batic heating or cooling. Since c 2 cx p^~ al , the pre-shock 
sound speed will in this case be a function of Se, 

Cs = C Si avg (8 I + 1) 1/3 , (7) 

for our 7 = 5/3 gas. In eq. ([7]), c s , avg is the "cosmic av- 
erage" sound speed (the sound speed of the intergalactic 
medium at a density equal to the cosmic mean at the 
epoch of interest), which we will again take to be equal 
to 15 km/s. 

In order to make further progress and be able to calcu- 
late measures of the statistical properties of the popula- 
tion shocks in this approximation, we need an analytical 
model for the environment of collapsed structures. For 
this purpose, we will use the double distribution (DD) 
of collapsed structures with respect to mass and local 
overdensity (Pavlidou & Fields 2005, hereafter PF05). 
In PF05, we defined the "local environment" of a col- 
lapsed structure through the clustering scale parameter, 
p. The clustering scale parameter is a free parameter in 
the double distribution model, and is defined so that the 
"environment" of an object of mass m be a surrounding 
region in space which encompasses mass (3m. Clearly, as 
{3 — ► 1 the environment of an object is restricted to be 
closer to the object itself, and for (3 — 1 it encompasses 
solely the collapsed object itself. Conversely, as (3 — > 00, 
the environment of the object extends further and fur- 
ther away from the object, to eventually encompass the 
whole universe. 
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With this definition of the "local environment", and 
using the same random walk formalism which yields the 
Press-Schechter mass function for collapsed objects, we 
found the DD to be 



dn 



PmS) S c (a) - St 



dmdSt m 



2tt 



dS_ 

dm 



exp 



(s c (a) - Stj 
2[S(m)- S(f3m)} 



exp 







2S(0m) 


— exp 



(5 e -2S c (a)) 2 
2S{f3m) 



[Stfm)} 1 / 2 [S(m) - S((3m)f /2 

In Eq. (|8|), S (to) is the variance of the linearly ex- 
trapolated overdensity field and St is the local linearly 
extrapolated overdensity (or underdensity) , which is re- 
lated to the true (calculated from the spherical evolution 

model) overdensity St of the environment sphere includ- 
ing the local object through the exact relations given in 
PF05 or through the useful approximation (accurate at 
a better than 2% level throughout its domain for all cos- 
mologies of interest) 



(8) 



D(a) 6c 



l-(l + &) 



-l/5 c 



(9) 



where D(a) is the linear growth factor in the cosmology 
of interest and S c is the linear overdensity threshold for 
collapse in the same cosmology. The quantity Si is in turn 
related to the overdensity of interest, St (the overdensity 
of the environment sphere excluding the central object), 
through 



St = 



Q3-i)(i + fr)(i + <y 

[3(1 + S c ) - (1 + S e ) 



- 1 



(10) 



Because in the problem of cosmic accretion shocks we 
are interested in the properties (density and sound speed) 
of the material right outside the shock surface, we will 
adopt a small value for the clustering scale parameter, 
= 1.1. Our results, however, are not sensitive to the 
exact value of [3 since, as we found in PF05, the prop- 
erties of the double distribution (when calculated as a 

function of St rather than St) depend only mildly on [3 
for small values of (3. 

3.3. Model 3 (full model): model 2 + effect of filaments 

In this third variation of our model, we will expand the 
doublc-distribution-based model 2 to include the effect 
of filament preheating and compression. Most virialized 
structures reside inside filaments, and hence they are ac- 
creting gas which has already compressed and preheated 
in filament shocks. A complete description of the ef- 
fect of filaments on the properties of cosmic accretion 
shocks would require an analytic model for the filamen- 
tary structure of the universe accounting for a distribu- 
tion of different filament temperatures and densities, as 
well as their evolution with redshift, a problem which we 
will address in detail in an upcoming publication. Here, 
we adopt a simplified prescription which can provide a 
first assessment of the degree to which filaments alter the 
energetics and evolution of cosmic accretion shocks. 



We will assume that all structures with environmen- 
tal overdensity St > reside inside filaments 3 . These 
structures, instead of accreting gas of mean sound speed 
Cs.avg, will be accreting filament-preheated gas of mean 
sound speed c Sj fii ~ 45km/s, corresponding to a (conser- 
vative) mean filament temperature Tgi « 10 5 K (e.g. Cen 
& Ostriker 1999, Dave et al. 2001; Fang et al. 2002). In 
addition, 1 + St will now represent the compression factor 
with respect to a mean filament density pai — Pb(l + ^fil) 
where (again, conservatively) 4 1 + <5ei ~ 10. Hence, in 
model 3 the pre-shock sound speed will be 



1 



,1-/3 



1 



,l./3 



St <0 
St > o 



c J C s,avg (St 
\ C a fil (Si - 

while the pre-shock density enhancement will be 

1 + Si , St<0 
(1 + S m )(l + St) , 5t>0 ■ 
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(11) 



(12) 



The distribution of number density of accretors with re- 
spect to mass and local overdensity St will again be as- 
sumed to be well described by the double distribution of 
collapsed structures. 

3.4. Quantitative measures of the properties of cosmic 
accretion shocks 

The specific properties of the population of cosmic ac- 
cretion shocks which we wish to investigate can be clas- 
sified into two categories. First of all, we are inter- 
ested to see how the shock-processed mass and energy 
is distributed among shocks of different Mach numbers, 
and how environmental effects alter this distribution. 
Such Mach number distributions are of particular im- 
portance to assess the role of cosmic accretion shocks as 
potential sites of particle acceleration, since the highest- 
energy particles can be accelerated only in the strongest 
(M ^> 1) shocks. The quantities we will calculate to 
address this question are: 

• The "mass current distribution" with respect to 
Mach number, dJ/d In A4, defined as the comov- 
ing mass current density crossing shock surfaces 
of logarithmic Mach number between ln.M and 
InM + dlnM, with units of M Q yr" 1 Mpc~ 3 . In 
the case of model 1, it is given by 



dJ 



dlnM 



M 



dn dm 



dm dM 



■h 



(13) 



3 This is a conservative assumption. In reality, even structures 
predicted by the double distribution model to reside inside under- 
densities may be found inside filaments, as neighboring voids may 
sufficiently compress their surrounding regions. 

4 Note that our prescription does not account for filament evo- 
lution with cosmic time. Of course, filaments do evolve with red- 
shift, and a more detailed calculation, based on a model account- 
ing for filament evolution as well as different filament overdensities 
and temperatures is desirable and currently pursued. However, 
the error introduced by not accounting for filament evolution is 
expected to be relatively small, as filaments primarily evolve by 
shock-heating and compressing increasing quantities of diffuse gas, 
rather than heating it and compressing it to significantly differ- 
ent densities and temperatures (see e.g. distributions in Dave et 
al. 200f, which change in amplitude rather than shift in location 
with changing redshift). This however is not expected to signifi- 
cantly affect the properties of gas accreted by virialized structures: 
filaments may include smaller quantities of gas at high redshifts, 
but structures are also smaller in size and most of them are still 
expected to be found inside filaments. 
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where m — m(M.,z) is given by eqs. (Q} and 
with c s i = 15km/s; dn/dm is the Press-Schecter 
mass function, given by eq. ©; and J\ is given by 
eq. ([3]) with S s = 0. In the case of models 2 and 3, 
it is given by 



d.J 



dlnM 



M 



d8 J ^ U ^ m 
,w=_x ^ 1 dmdde dM. 



(14) 



5^ ,2 



where dn/(dmdSi) is the double distribution of 
cosmic structures, given by eq. (|8|); m(M,S, z) is 
again given by eqs. fl} and @, with c s i given by 
eqs. §7$ and (fTTj) for models 2 and 3 respectively; 
and S s = 8 1 for model 2, while (5 S given by eq. (fT2"j) 
for model 3. 

• The "kinetic power distribution" with respect to 
Mach number, dP/d\nM, defined as the comov- 
ing kinetic power density crossing shock surfaces 
of logarithmic Mach number between lnA4 and 
InM + dliiM, with units ergs -1 Mpc -3 . It is 
given by 



dP 



dhiM 
for model 1 and by 
dP 



M 



dn dm 



dm dM 



Pi 



d\nM 



M 



dn dm 



s t =-i amade dM 



(15) 



(16) 



for models 2 and 3. In both cases, P\ (M , z) is given 
by Eq. ©. 

Second, we would like to follow the cosmic evolution of 
the energetics of cosmic accretion shocks, tracing their 
mass and energy processing and their impact on the 
intergalactic medium. In particular, we wish to calcu- 
late the energy amount processed by all cosmic accretion 
shocks per unit time and by a certain redshift, as well as 
the degree to which cosmic baryons have been affected by 
cosmic accretion shocks by some specific cosmic epoch. 
The quantities we will use to quantify these questions 
are: 

• The "integrated mass current density" , J, which is 
the comoving mass current density crossing shock 
surfaces of any Mach number at a given cosmic 
epoch, with units of M Q yr _1 Mpc -3 . 

• The "integrated kinetic power density" , P, which is 
the comoving kinetic power density crossing shock 
surfaces of any Mach number at a given cosmic 
epoch, with units of ergs -1 Mpc~ 3 . 

• The "cumulative processed kinetic energy density" , 
J t Pdt, which is the total kinetic energy density 
processed by shocks of any Mach number since 
some initial cosmic epoch t- u with units of eV per 
baryon in the universe. 

The second set of quantities will be calculated as a 
function of redshift, and for each of the models 1, 2 and 3, 
in order to assess the effect of each environmental factor 
on their behavior. 



4. RESULTS 

Throughout this section we will assume a flat, 
vacuum+matter universe with Wilkinson Microwave 
Anisotropy Probe (WMAP) parameters a s — 0.84, h = 
0.71, n m = 0.27 and Sl h = 0.04 (Spergel et al. 2003). 
We focus on post-reionization redshifts, hence we adopt 
a cosmic average sound speed of 15km/s, corresponding 
to a temperature of ~ 10 4 K, for a fully ionized plasma 
with fi — 0.59 (25% He by mass). 

4.1. Mach number distributions and the effect of the 
environment 

The distributions with respect to Mach number of mass 
current and kinetic power are plotted in figs. [T] and [2] re- 
spectively. In both figures, the solid line corresponds 
to model 1 (Press-Schechter-based, no environmental ef- 
fects); the dashed line corresponds to model 2 (double- 
distribution-based, includes the effect of the primordial 
density field but not the effect of filaments); and the 
dot-dashed line corresponds to model 3 (includes effects 
of both primordial density field as well as effects of fila- 
ments). The left panel shows the results for z = 3, while 
the right panel corresponds to the present-day universe. 

The results of the Press-Schccter-based, no- 
environmental-effects model can be interpreted in 
a straight-forward way. In this case, accretor mass and 
Mach number have a one-to-one correspondence, with 
larger accretor masses resulting to higher— Mach-number 
accretion shocks. As cosmic time progresses, the mass of 
the largest collapsed objects increases, as more and more 
massive structures have had enough time to collapse. 
Hence, the strongest accretion shocks become stronger 
with decreasing redshift, and the distributions of mass 
current and kinetic energy move to higher Mach numbers 
with increasing time. The kinetic power distribution is 
more strongly suppressed at low Mach numbers (and 
therefore low accretor masses) than the mass current 
distribution. This result can be immediately verified 
by simple analytic arguments. For masses high enough 
that the primordial density fluctuation power spectrum 
can be regarded as a power law but low enough that 
the exponential mass cutoff is not affecting the results, 
dn/dm oc m~ 2 . In addition, M oc m 1 / 3 (in the high-A'J 
limit), while r v tx m 1 / 3 . Hence, J\ oc m and P\ oc m 5 / 3 . 
Equations (fT3|) and (fT5]) then give for the low-mass 
dependence of the mass current and kinetic power 
distributions, 



dJ 



d\nM 
dP 



ocm° oc M° 
ocm 2//3 oc M 2 



d^M- ' (17) 
accounting for the difference in the low-mass behavior of 
different distributions. Note that the low-mass suppres- 
sion is somewhat stronger than what predicted by the 
simple arguments above, due to the deviation of A4(m) 
from m 1 / 3 for M. — > 1. 

It should be pointed out that the result that it is the 
highest-Mach-number shocks that process most energy 
refers to the energy processed by accretion shocks alone 
(which are the shocks treated in this work). Merger 
shocks are generally expected to dominate the overall en- 
ergetics of cosmic shocks (e.g. Miniati et al. 2000; Gabici 
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Fig. 1. — Mass current distribution (spatial density of mass accretion rate per logarithmic Mach number interval, dj/d\n A4)for z = 3 
(left panel) and 2 = (right panel). The units of the vertical axes are Mq yr — 1 comoving Mpc - 3 . Solid line: model 1 (no environmental 
effects); dashed line: model 2 (effect of primordial density fluctuations); dot-dashed line: model 3 (effects of primordial density fluctuations 
and filaments). 
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Fig. 2. — Kinetic power distribution (spatial density of kinetic power processed by accretion shocks per logarithmic Mach number interval, 
dP/d\n A4)for 2 = 3 (left panel) and 2 = (right panel). The units of the vertical axes are g s" 1 comovingMpc -3 . Solid line: model 1 
(no environmental effects); dashed line: model 2 (effect of primordial density fluctuations); dot-dashed line: model 3 (effects of primordial 
density fluctuations and filaments). 



& Blasi 2003; Ryu et al. 2003), because the density of the 
(already virialized) gas they process is much higher (by 
a factor between ~ 20 and 200 on average, depending on 
the environment in which accretion shocks lie) than the 
density of the gas processed by accretion shocks. How- 
ever, because the virialized gas relevant to merger shocks 
has also been preheated to much higher temperatures 
compared to the gas accreted by accretion shocks, the 
Mach numbers associated with merger shocks are signif- 
icantly lower. Hence, accretion shocks are expected to 
dominate the strong shock energetics. 
The effect of taking into account the environmental 



density fluctuations in the primordial density field using 
the double distribution of collapsed structures is demon- 
strated by the curves corresponding to model 2, and is 
two-fold. On the one hand, every mass bin is spread 
out to a larger M range, as each accretor mass now cor- 
responds to a distribution of Mach numbers, reflecting 
the associated distribution of environmental overdensi- 
ties. On the other hand, the local density of the material 
immediately outside the accretor is now correlated with 
the accretor mass. Higher-mass objects, which have a 
dominant contribution to the kinetic power and, to a 
lesser extent, the mass processed by accretion shocks, 
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tend to reside in higher-mass environments. As a result, 
the overall amplitude of both the kinetic power and the 
mass current distribution is increased. 

A most striking environmental effect in the case of the 
kinetic power distribution is that, as redshift decreases, 
a second peak separates out in model 2, which is absent 
in the base model. The presence of this second, high-Al 
peak is due to the gradual shift of increasingly massive 
structures towards underdense environments. The lower- 
A4 peak is, conversely, due to the higher-mass structures 
which are embedded inside higher-overdensity and hence 
higher-temperature environments. 

The additional effect of the filamentary structure of 
the universe is demonstrated by the curves correspond- 
ing to model 3. The inclusion of filaments results to a 
shifting of both the mass current and the kinetic power 
distributions towards lower Mach numbers (due to the 
preheating of the pre-shock gas by filament shocks) and 
higher amplitude (due to the additional density enhance- 
ment of pre-shock material inside filaments). At low red- 
shifts, the amplitude of both distributions shifts by an or- 
der of magnitude, reflecting our assumed mean filament 
overdensity. Although our recipe for the inclusion of the 
effect of filaments in model 3 is only a first approxima- 
tion to the analytical treatment of the issue, there are 
two robust conclusions we already draw from this anal- 
ysis. First, the effect of filaments is at least of the mag- 
nitude derived here, since all of our assumptions were 
chosen so that they err on the side of underestimating 
both the filament-associated preheating as well as den- 
sity enhancement. Hence, our results emphasize the need 
for a detailed model for filament-associated processes and 
their cosmic evolution, which have been shown to be an 
important factor in determining the statistical proper- 
ties of cosmic accretion shocks. Second, the direction of 
the changes in the mass current and kinetic power dis- 
tributions due to filaments will be towards lower Mach 
numbers and higher amplitudes, if anything at a higher 
degree than predicted here. Even with more sophisti- 
cated models, this result is not expected to change qual- 
itatively. 

Note that the mean filament temperature we have used 
(T ~ 10 5 K) is rather on the low-end of the temperature 
distribution of filaments in the present-day universe as 
found in numerical simulations. Since we do not account 
for redshift evolution of the mean filament temperature, 
we have chosen this low value to ensure that we do not 
overestimate the effect of filament preheating on accre- 
tion shocks at higher redshifts. A more realistic mean 
filament temperature of T ~ 10 6 K at low redshifts (right 
panels of Figs. [T] and [J) would roughly result to a hor- 
izontal translation of the model 3 curve towards lower 
Mach numbers by a factor of ~ 3. 

4.2. Integrated quantities and evolution with cosmic 

epoch 

Figure [3] shows the evolution of the integrated mass 
current J (left panel) and integrated kinetic power P 
(right panel) for shocks of any Mach number, and for 
redshifts between 6 and 0. The solid line corresponds 
to model 1, the dashed line to model 2, and the dot- 
dashed line to the "full" model 3. The thick solid line 
corresponds to the supernova energy input to the inter- 
galactic medium (which is discussed in detail in the next 



section). 

The integrated mass current peaks at z ~ 2 in the case 
of model 3, while the models 1 and 2 mass current histo- 
ries peak closer to z ~ 3. The reason for this difference is 
that the filament preheating and the associated decrease 
of accretion shock Mach number for given accretor mass 
result to a certain loss of power from the mass current 
distribution, as the lowest mass objects can no longer 
satisfy A4 > 1. The effect is more pronounced at higher 
redshifts, when the largest collapsed objects are less mas- 
sive than the largest objects in the local universe. As a 
result, the evolution of the mass current is more severe 
at high redshifts for model 3, and the integrated mass 
current peaks at a lower redshift than in the other two 
models. 

The same effect is present but less pronounced in the 
case of the integrated kinetic power: the high-z tail of 
model 3 falls more sharply with increasing redshift than 
the tail of models 1 and 2, however the effect is not strong 
enough to affect the location of the peak, which is found 
at z <~ 1 for all three models. 

The difference in the location of the peak in redshift be- 
tween integrated mass current and kinetic power is due to 
the difference in the mass dependence of the associated 
distributions. The kinetic power distribution depends 
on mass much more strongly than the mass current, 
and thus the integrated kinetic power history "prefers" 
lower redshifts, where the largest collapsed structures 
(which process most of the kinetic power crossing ac- 
cretion shocks) are more massive. For the same reason, 
the "loss of power" due to the presence of filaments is 
much milder in the kinetic power history: the lower-mass 
structures, which are in danger of "dropping out" from 
the population of structures hosting shocks in a higher- 
temperature environment, do not contribute much to the 
kinetic power distribution to begin with. 

At redshifts close to their respective peaks, the in- 
tegrated mass current and kinetic power as calculated 
by model 2 are enhanced by about an order of magni- 
tude with respect to the no-environmental-effects model 
1. Similarly, the predictions of model 3 (filaments + pri- 
mordial overdensities) are about an order of magnitude 
enhanced with respect to those of model 2 (primordial 
overdensities only). 

4.3. Comparison of accretion shocks with other 
intergalactic medium energy inputs 

In this section, we will discuss how the kinetic en- 
ergy processed by cosmic accretion shocks compares with 
other energy inputs and characteristic energy scales of 
the intergalactic medium. 

In the right panel of fig. Owe have plotted, along with 
the integrated kinetic power of cosmic accretion shocks 
as predicted by models 1, 2 and 3, the energy input of 
type II supernovae as a function of redshift (thick solid 
line) . The supernova energy history was derived from the 
Cole at al (2001) fit to their dust-corrected measurements 
of the cosmic star formation rate, assuming a Salpeter 
mass function with a supernova progenitor mass cutoff 
of 8Mq, and adopting a supernova explosion energy of 
10 51 erg. The accretion-shock-processed energy is higher 
than the energy output of supernovae for all redshifts 
< 3. In addition, at low redshifts the star formation rate 
decreases with decreasing redshift much more strongly 
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Fig. 3. — Integrated mass current J (left panel) and integrated kinetic power P (right panel) as a function of redshift. Solid line: model 1 
(no environmental effects); dashed line: model 2 (effect of primordial density fluctuations); dot-dashed line: model 3 (effects of primordial 
density fluctuations and filaments). The thick solid line in the right panel represents the energy output of supernovae as derived from the 
Cole et al. (2001) fit for the cosmic star formation rate. 



than the accretion shock processed power, and as a result 
in the local universe the energy input to the intergalactic 
medium due to accretion shocks is expected to be (at 
least) one order of magnitude higher than the energy 
output of supernovae. 

^From the point of view of particle acceleration, de- 
spite the fact that the efficiency with which accretion 
shocks may accelerate high-energy particles is largely un- 
constrained, the energy available for particle acceleration 
is much larger than in the case of supernovae, and in 
this respect the potential of cosmic accretion shocks as 
sites of cosmic ray acceleration is once again seen to be 
promising. 

Finally, in Fig. [4] we plot the cosmic history of J Pdt 
(the cumulative processed kinetic energy) in units of eV 
per baryon in the universe (as opposed to per shocked 
baryon). Again, the solid line corresponds to model 1, 
the dashed line to model 2 and the dot-dashed line to 
the "full" model 3. The horizontal line in this plot corre- 
sponds to 13.6 eV per baryon. ^From the location of the 
intersection of the horizontal line with the J Pdt curve, 
we can conclude that by redshift z ~ 3.5 (for model 3), 
the energy processed by accretion shocks alone would 
have become comparable to the energy needed to reion- 
ize the universe 5 . 

5. COMPARISON TO PAST WORK 

5.1. Comparison to baryon energetics from cosmological 
simulations 

If we integrate J (the mass current through shocks of 
any Mach number) as predicted by model 3 (dot-dashed 
line in left panel of Fig. [3|) over time, we find that by 

J This is an order-of-magnitude estimate, meant to give a feeling 
of the amount of energy processed by shocks as compared to other 
energy inputs in the ION!. If one wanted to consider shocks as an 
actual reionization mechanism, a detailed modeling of the reion- 
ization process would be required, see e.g. Miniati et al. 2004. 




z 



Fig. 4. — Cumulative kinetic power J Pdt (in units of eV per 
baryon in the universe), as a function of redshift. Solid line: model 
1: dashed line: model 2; dot-dashed line: model 3. The horizontal 
dotted line corresponds to 13.6 eV per baryon. 



the present cosmic epoch, 51% of the baryons in the uni- 
verse have been processed by accretion shocks, and are 
therefore residing within structures of virial temperature 
roughly > 10 4 K (mass > 10 10 Mq). We can compare this 
value to the findings of cosmological simulations. If we 
assume that 1 /2 of the baryons of temperatures between 
10 4 — 10 7 K belong to collapsed objects of correspond- 
ing virial temperatures, and we add these baryons to 
the "condensed" (cooled and belonging to galaxies) and 
"hot" (T > 10 7 K and belonging to high-mass collapsed 
objects) baryons, we get a total baryon fraction which 
is between 35 — 60%, depending on the specifics of each 
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simulation (see, e.g. simulations in Dave et at 2001). 

As far as the cumulative energy processed by accretion 
shocks J Pdt is concerned, the prediction of model 3 for 
by z = is P ks 2 keV per baryon in the universe (dot- 
dashed line in Fig. 0]) . Note that this is not a prediction 
for the mean per baryon in the universe energy today, 
since our calculation does not account for any cooling 
losses or any other energy inputs to the IGM. Still, it is 
interesting to note that it is within a factor of ~ 2 of 
the value found in cosmological simulations by e.g. Cen 
& Ostriker (1999), who (accounting for cooling as well 
as star formation feedback) calculate a density-weighted 
average temperature (the quantity representing the mean 
per-baryon-in-the-universe energy) of ~ 10 7 K at z = 0. 

It is therefore very encouraging to see that, given the 
simplifying assumptions and idealizations of our model 
(spherically symmetric infall, Press-Schechter mass func- 
tion, absence of time evolution in filament temperature 
and density), our results are in reasonable agreement 
with the findings of cosmological simulations, indicat- 
ing that our most sophisticated model (model 3) has 
captured the essential elements of the gravity-driven, 
diffuse-matter-accretion-associated energetics of cosmic 
baryons. More importantly, we see that our other two 
models, which differ 1-2 orders of magnitude in their pre- 
dictions for the cumulative quantities discussed above, 
are not sufficient to reproduce the results of cosmological 
simulations. This fact is suggestive of the importance 
of filaments not only in processing baryons in their own 
filament shocks and giving rise to the warm-hot inter- 
galactic medium, but also in modifying the energetics of 
accretion shocks. 

5.2. Comparison to simulations of cosmic shocks 

Figure [5] shows two quantities commonly used in the 
cosmic shock literature to describe the properties of pop- 
ulations of different types of cosmic shocks. In the left 
panel, we plot the distribution of the comoving num- 
ber density of accreting structures per logarithmic Mach 
number interval of their respective accretion shocks, 
dn/dhiA4, with units number of structures per comov- 
ing Mpc 3 , for the three models we have presented in this 
work and for z = 0. In the right panel, we plot the distri- 
bution of the comoving shock surface area per logarith- 
mic Mach number interval per volume under considera- 
tion, dS/dhxM. with units of Mpc -1 (ratio of shock sur- 
face over space volume) . Only structures with virial tem- 
peratures > 10 4 K (masses > 8 x 10 9 M©) are plotted. As 
expected, these distributions are dominated by the low- 
mass, low-Mach-number structures (the same analytic 
scalings used in §4.11 give, for the low-mass behavor of 
our no-environmental-effects model 1 dn/dhiA4 oc mT x 
and dS/d\nAi oc m -1 / 3 ). The number distribution 
in the Press-Schechter-based model (model 1) mono- 
tonically increases for decreasing Ai since in this case 
there is a one-to-one correspondence between mass and 
Mach number. Hence the number distribution of objects 
simply follows the Press-Schechetr mass function modu- 
lated by the mass-Mach conversion. In the case of the 
double-distribution-based models (models 2 and 3), the 
mass cutoff, in combination with the distribution of pre- 
shock sound speeds for each accretor mass described by 
the double distribution, results in a number distribution 



which peaks at M < 4, a position which is defined by 
the combination of the virial temperature at the mass 
cutoff (which is the most populated available mass bin) 
and the sound speed corresponding to the most proba- 
ble environment at that particular mass. Note that the 
difference in both the number and surface distributions 
between models 2 and 3 (without and with the effect of 
filaments) is very small. This is because we have used the 
double distribution of cosmic structures to predict which 
structures are located inside filaments, and the double 
distribution predicts that most of the low-mass cosmic 
structures (the ones dominating the number and surface 
distributions) are found inside underdense regions, and 
hence not inside filaments. 

For the case of accretion shocks which we have studied 
here, we can then draw the following conclusions: 

1. The energetics of these shocks are determined by 
the few, high-Mach-number, high-mass structures, 
and are not significantly affected by the majority 
of structures, which have low masses and low ac- 
cretion shock Mach numbers. 

2. Our results are in natural agreement with the find- 
ings of Ryu et al. (2003), who performed a series of 
simulations with increasing mass and spatial res- 
olution, and found convergence for their results 
on kinetic power distribution (dominated by high- 
mass structures) but not for the surface distribu- 
tion (dominated by low-mass structures). However, 
since the physical impact of shocks to their envi- 
ronment is better represented by the kinetic power 
distribution, which is not sensitive to the behavior 
of low-mass structures, the lack of convergence in 
the number-dominated distributions in cosmologi- 
cal simulations should not decrease confidence in 
the relevance and robustness of the conclusions of 
such cosmological simulations relating to the effect 
of accretion shock on the intergalactic medium. 

Finally, we should caution the reader that a direct com- 
parison between our Fig. ([5]) and similar plots in cos- 
mic shock literature is not straight-forward - the shocks 
treated here are exclusively accretion shocks, and they 
do not include merger or filament shocks. When shock 
surfaces in the universe are labeled as internal or exter- 
nal, accretion shocks are split between the two types: 
internal shocks can be merger shocks, or they can be 
accretion shocks associated with structures found inside 
filaments; external shocks on the other hand can be ac- 
cretion shocks associated with structures which are not 
located inside filaments, or they can be filament shocks. 

5.3. Comparison to other analytic models of cosmic 
shock populations 

Press-Schechter extensions have been used extensively 
to model different "families" of large-scale cosmic shocks, 
with emphasis primarily on the emission properties of 
such populations, rather than the properties of the 
shocks themselves. 

Waxman & Loeb (2000) assumed a Press-Schechter 
mass function to describe an underlying population of ac- 
cretors hosting particle accelerating shocks which would 
produce, through their radio and gamma emission, fluc- 
tuations in the radio and gamma-ray background. All 
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Fig. 5. — Left panel: distribution of the comoving number density of accreting structures per logarithmic Mach number interval of their 
respective accretion shocks, dn/dlnAi, with units number of structures per comoving Mpc . Right panel: distribution of the comoving 
shock surface area per logarithmic Mach number interval per volume under consideration, dS/dln M with units of Mpc -1 . Solid line: 
model 1; dashed line: model 2; dot-dashed line: model 3. All curves correspond to z = 0. Only structures with masses > 8 X 10 9 Mq are 
plotted. 



shock quantities required for their emissivity and fluctu- 
ation calculations (mass accretion rate, shock tempera- 
ture, shock radius) were derived as functions of the accre- 
tor velocity dispersion up to dimensionless constants, left 
to be calibrated against simulations. Of all the analytic 
shock models discussed in this section, the Waxman & 
Loeb (2000) model is the one with the most similarities 
with the "base" model we use in this work (the model 
containing no environmental effects). The two primary 
similarities are that a) we also take the underlying ac- 
cretor distribution to be described by a Press-Schechter 
mass function and b) the mass-scaling of the shock prop- 
erties is the same over a large range of masses. In our 
case, the shock properties are obtained the tempera- 
ture jump across the shock, which involves the mass- 
dependent virial temperature of each object and some 
environmental temperature. At the large-temperature- 
jump/high-Mach- number limit, the mass dependencies 
of shock properties in our model agree with those of 
Waxman & Loeb (2000), as the velocity dispersion and 
the free-fall velocity around an object of a certain mass 
(which is encoded in the accretion shock Mach number 
provided that Mach ^> 1, or, equivalently, that the gas 
pressure effects in the pre-shock gas are negligible) nec- 
essarily have the same mass dependence. 

However, our "base model" describes different aspects 
of the accretion shock population that Waxman & Loeb 
(2000), as they calculate the shock gamma and radio 
emission properties while we calculate the mass current 
and shock energetics, exhibit their dependence on Mach 
number, calculate the cosmic distribution of mass and 
energy flows as a function of Mach number, and explicitly 
construct histories for the accreted mass and energy. 

Fujita & Sarazin (2001) used a version of the Lacey 
& Cole (1993) merger formalism (a Press-Schechter ex- 
tension) to describe mergers and accretion in a unified 



way (in their picture, accretion is the sum of all mi- 
nor mergers for which the relative mass increment of the 
accreting halo does not exceed some threshold). They 
as well use their model to obtain the nonthermal emis- 
sion properties of merging and accreting clusters, and 
they compare how such properties differ between clusters 
which have recently undergone some major merger and 
clusters which have been only accreting through minor 
mergers. Treating accretion as a subcase of the general 
merger picture associated with hierarchical structure for- 
mation is a widely used technique. It was adopted also by 
Gabici & Blasi (2003a, b) who also used a semi-analytic 
model based on the Lacey & Cole (1993) merger formal- 
ism to calculate the Mach numbers of shocks resulting 
from mergers, the particle acceleration properties of such 
shocks, and the contribution of merger-shock-accelerated 
particles to nonthermal radiation from large-scale struc- 
ture (including cluster radio halos and a contribution to 
the gamma-ray background) 6 . 

Such a mergers-only approach is fundamentally rather 
than just formally different than the one we use here. 
First of all, since in the mergers-only picture all matter 
is locked up in collapsed objects of various masses, the 
pre-shock material has a mixture of temperatures, char- 
acteristic of the virial temperatures (and hence masses) 
of the accumulated halos, rather than of the diffuse en- 
vironment of the object. In addition, the mass accretion 
rate for a halo of mass m in the mergers-only picture is 
not oc m as in the simple, gravitationally-driven accre- 
tion of diffuse gas picture we adopt, but is dependent on 
the index n of the power spectrum of initial perturbations 

6 Note however that Gabici & Blasi (2004) did make an explicit 
distinction between mergers and diffuse gas accretion when assess- 
ing the detcctability of gamma rays from galaxy clusters; in this 
work, they considered all accretion shocks to have JVi 2> 1, while 
the scaling of shock energy with accretor mass was similar to that 
of Waxman & Loeb (2000). 
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(and is oc m 7 / 6 for n = 1). This difference stems from the 
different physical processes followed in each case. In the 
mergers-only, Press-Schechter-extension picture, growth 
of structure, whether through accretion or through merg- 
ers, occurs due to the gradual turnaround and collapse 
of increasingly large scales. This process is caused by the 
gravitational enhancement of perturbations already im- 
printed in the primordial density field. The gravitational 
effect of a collapsed object within such a larger collaps- 
ing structure is simply to contribute to the mean den- 
sity of the collapsing region. In our picture, accretion is 
driven locally, through the gravitational attraction that 
a collapsed structure exerts on its surrounding diffuse 
gas. Each collapsed structure would accrete gas even if 
it lived inside an otherwise homogeneous universe. In 
nature, both processes contribute to the growth of any 
particular collapsed structure, and it is hence important 
to construct analytic models based on both pictures, and 
gauge the relevant importance of each process through 
comparison with observations and simulations. 

More recently, Inoue & Nagashima (2005) used the 
Somerville & Kolatt (1999) multiple merger model to fol- 
low the population of cosmic shocks, calculate the associ- 
ated gamma-ray emission and investigate the possibility 
to use gamma-ray observations as a tool for the study 
of the warm-hot intergalactic medium. The Somerville 
& Kolatt (1999) merger model is an improved version of 
a Press-Schechter-based semi-analytic merger-tree con- 
struction algorithm. Its major advantage with respect 
to other algorithms for construction of merger trees of 
present-day objects is that it predicts a distribution of 
mass progenitors which is consistent with the extended 
Press-Schechter predictions while it enforces strict mass 
conservation. This model accounts for some amount of 
diffuse accreted matter, which is however calculated not 
by gravitational arguments, but again as that amount 
of matter corresponding to progenitors less massive than 
some pre-set cutoff. 

Finally, Furnaletto & Loeb (2004) used a modified 
Press-Schechter formalism to describe large-scale shocks 
that may appear when overdense perturbations reach and 
exceed their turnaround point. The shocks they con- 
sider are located at the boundaries between regions which 
have reached their turnaround point, and the background 
(still expanding) universe. The shock velocity is derived 
from a scaling argument to be approximately equal to 
the size the shocked region would have had if it had 
expanded with the Hubble flow, divided by the age of 
the universe. All shocks in this work are assumed to 
be strong, and Mach-number dependencies are not in- 
vestigated. ^From a physical standpoint, the shocks de- 
scribed by these authors occur much earlier in the evo- 
lution of each structure than the accretion shocks we 
describe here, and they are the shocks responsible for 
heating up the bulk baryonic component of collapsed ob- 
jects from mean intergalactic medium temperatures to 
the virial temperature of structures. In our shock clas- 
sification scheme, the Furlanetto & Loeb (2004) shocks 
would be the gravitationally-driven component of fila- 
ment shocks. 

To summarize, our model of accretion shocks is dif- 
ferent from and complementary to all of the exiting 
treatments. Our focus on the role of diffuse gas accre- 
tion is qualitatively and quantitatively different from the 



mergers-only approaches, and these differences can be 
exploited to test against the physical reality in which 
both accretion modes occur. Our model is most sim- 
ilar to that of Waxman & Loeb (2000), but we extend 
the Press-Schechter treatment to include both the effects 
of environment, and to address the effects of preheating 
by filament shocks. Finally, our work differs from much 
of the existing literature by focusing on the shock en- 
ergetics, and to our knowledge is the first to explicitly 
calculate the net cosmic processing of both energy and 
mass through accretion shocks over cosmic history. 

6. DISCUSSION 

In this paper, we have investigated analytically the 
properties of cosmic accretion shocks around collapsed 
structures. We have calculated the mass and kinetic en- 
ergy currents as a function of Mach number for accretion 
shocks on individual structures of arbitrary mass at any 
epoch. Using this, we have computed the cosmic distri- 
bution of mass and energy flows as a function of Mach 
number, and we have calculated the evolution of these 
quantities as a function of rcdshift. 

We have found that environmental factors play a major 
role in shaping the qualitative and quantitative proper- 
ties of accretion shocks. We demonstrated the impact 
of environmental effects by exploring three models for 
the cosmic shock population. The first used the Press- 
Schechter mass function to describe the underlying pop- 
ulation of collapsed, accreting objects. All such objects 
were assumed to accrete material of the same density 
and temperature. This was our "control" model, which 
did not include any environmental effects. The second 
model used the double distribution of collapsed struc- 
tures (Pavlidou & Fields 2005) to describe the distri- 
bution of accreting objects with respect to both their 
mass and local environment overdensity or underdensity. 
The overall mass distribution of objects is the same as 
in the first model, as the double distribution integrates 
to the Press-Schechter mass function. This model then 
only included environmental effects caused by fluctua- 
tions in the primordial density field. Our third model 
was a refinement of model 2, which considered all ob- 
jects predicted by the double distribution to reside inside 
overdensities to accrete gas preheated and compressed in 
cosmic filaments. 

The inclusion of primordial density fluctuations re- 
sulted to broader distributions with respect to Mach 
number of the mass current and kinetic power processed 
by accretion shocks, and an overall increase in the am- 
plitude of such distributions, owing to an increase in the 
density of accreted material by the most massive objects, 
which dominate especially the processed kinetic power. 
Perhaps the most striking effect of accounting for this 
environmental factor is the separation of a second, high 
Mach number peak in the kinetic power distribution at 
low redshifts, due to an increasing number of higher- mass 
structures concentrating inside underdensities. 

The inclusion of filament preheating and compression 
further increased the amplitude of the mass current and 
kinetic power distributions, and shifted the distributions 
towards lower Mach numbers. A more refined filament 
model, the need of which is underlined by our results, is 
expected to enhance these trends even more. Our treat- 
ment of both environmental factors has been systemati- 
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cally conservative in estimating the impact of the effects 
on the shock properties. Consequently, we expect that 
these differences will be evident at least qualitatively in 
full numerical models. 

The integrated kinetic power processed by shocks 
peaks at a redshift of ~ 1, while the integrated mass cur- 
rent peaks at higher redshifts, z ~ 2 (for model 3). The 
effect of the local environment is to increase the overall 
level of the processed energy at peak redshifts by a factor 
of ~ 10 for each environmental factor included. 

Comparing the energy processed by cosmic accre- 
tion shocks to other energy inputs to the intergalac- 
tic medium, we find that the energy input of accretion 
shocks (as predicted by model 3) is higher that that of 
Type II supernovae for all z < 3, and it becomes more 
than an order of magnitude higher in the local universe. 
In addition, we found that energy processed by accretion 
shocks alone becomes comparable to the energy needed 
to reionize the universe by z ~ 3.5. 

Numerical and observational tests of our model will 
quantify the importance of accretion shocks and their 
prominence among cosmic shock mechanisms. Our em- 
phasis on accretion due to diffuse matter is complemen- 
tary to minor mergers, thus comparison with simulations 
will shed light on the relative importance of these two 
components in the ACDM scenario. To compare the 
present results with simulation is however neither trivial 
nor immediate, because simulations include merger and 
filament shocks as well, but moreover because accretion 
shocks as we have defined them do not uniquely map onto 
the classification schemes used in the existing literature. 
A detailed, self-consistent comparison with simulations 
will appear in future work. 

Of course, the question of the nature of cosmic shocks 



and the importance of accretion processes ultimately will 
be resolved observationally. Galaxy clusters provide a 
promising site to study accretion shocks (e.g. Lieu et 
al. 1996; Ensslin et al. 1998; Fusco-Femiano et al. 1999; 
Fujita & Sarazin 2001; Govoni et al. 2001; Miniati et 
al. 2001a,b; Bagchi et al. 2002; Berrington & Dermer 
2003; Feretti et al. 2004; Gabici & Blasi 2004; Keshet et 
al. 2004; Kocsis et al. 2005) As with simulations, there is 
a need to separate merger versus accretion components; 
this might be done spatially (e.g. Inoue et al. 2005), 
with accretion shocks dominating peripheral emission, 
particularly that due to inverse Compton. In addition, 
detection of 7— ray emission unambiguously associated 
with structure formation shocks, either by TeV (such as 
CANGAROO, HESS, MAGIC and VERITAS) or GeV 
(such as AGILE and GLAST) gamma-ray telescopes, will 
provide invaluable insight not only in the nature of the 
accretion process itself, but also the currently elusive sub- 
jects of the acceleration efficiency associated with cosmic 
shocks, as well as the properties of large-scale magnetic 
fields. 
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